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The Hexagon Problem 


The seven hexagons shown in the Figure are to be 
rearranged in position and rotation so that the pairs of 
numbers at the adjacent edges have no factor in common 
except 1. In the pattern as shown, this is true for the 
edges between hexagons 2 and 7, and 5 and 6, but is not true 
for the edges between hexagons 2 and 3, or 4 and 5. 


There are 61-60 possible arrangements of the hexagons 
(1.e., 33,592,320 arrangements) but not all of them need to 
be tested. For any given arrangement, there are 12 pairs 
of numbers to test for a greatest common divisor of 1, but 
again, not all 12 tests need to be made, since any failure 
vitiates the need for further testing on that arrangement. 
A flowchart for an efficient scheme for solving the hexagon 
puzzle will be given in a later issue. 


The problem is due to Lee Morgenstern; it has a 
unique solution. 


Beginning with this issue, problems presented for 

computer solution are numbered. The 8 problems so far 

presented are as follows: 

Problem 1. The 3X+1 Problem, in PCl-1 and PCl-6. 

Problem 2. Multiples of 3 Problem, in PC1l-5 and PC4-8, 

Problem 3. The Four 4's Problem, in PC2-3, 

Problem 4. The Zigzag Problem, in PC2-13. 
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Problem 5. The Fourway Problem, in PC3-1. 

Problem 6, Numbering the Fractions, in PC3-9. 
Problem 7. The AP Problem, in PC3+13. 

Problem 8. Cycle Length of Reciprocals, in PC4-13, 
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Log 5 
In 5 


Ss 


0.6989700043 360188047 8626110527 550697 32318101185379 ~ 


1.6094379124 341003746007 593332261876395256013542685 
2 .2360679774997896964.0917 36687 312762354406183596115 
1.709975946676696989 35310887254. 38601098680551105431 
1.3797296614612148 3239006 34642160176928556498779776 
1.2584989506418267 349927871711777138948186804815067 
1.17461894308801900591446 366569 189890167 37230187574 
1.016224591267 325635816890679158819446959800747 6627 
148 .41315910257660342111558004055227962 3487667 59388 
306 .01968478528145326274.131004 34 35606480 30070662807 
1.37 34007669450158608612719264449611486509995958997 


7888609052210118054117 285652827 8622967 3206435109023 
0047702789 306640625 


93326 3618503218878990089544723817169617091446371708 
0246217143397959669 109757756 34454440327 097881102359 
5949899 3032424. 2624215487521 35403239484.1520817203930 
7562344.1066613832515027 3995075985901831511100490796 
265113118240512514795933790805178271125415103810698 
3788544264811194698142286609592220176629104427984.56 
169448887 147466528006 3283684 52647429261829862165202 
79319528949 360711785066 36687410654 398055307 18136320 
59984482604 1954101213229629869502194514609904214608 
668 3612447929520 348268646176579269160474200659 36 389 
0417 37895822118 36507804555662844427 3925387517127854 
796781556 3464.037 1487768176689985539206 3601631001004 
21197 3674701749862626690747 2967 6253580 3929376233833 
98104692787455860525 3696441650390625 
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Book Reew 


A GUIDE TO FORTRAN IV PROGRAMMING (Second Edition) by 
Daniel D. McCracken, John Wiley & Sons, 1972, soft cover, 


288 pages, $6.95. 
Reviewed by Edward A. Ryan, Woodland Hills, California 


The catalog of the library of books that have earned 
the accolade of "classic" in their field would be a small 
one, indeed. McCracken's A GUIDE TO FORTRAN IV PROGRAMMING 
would most certainly be a member of that library. In 
the short span of time since its introduction, McCracken's 
book and its FORTRAN II predecessor have outsold all 
other texts on the subject of FORTRAN combined. More 
importantly, the book appears to have withstood the test 
of time, as is clearly evident in the recently published 
second edition. 


Anyone familiar with the first edition will immediately 
recognize the clear, crisp, readable style of writing that 
distinguishes this book from the stodgy, language manual 
style of so many of its competitors. The new edition is 
nearly double the size of the previous version. The 
organization and presentation of the material follows 
that of the earlier version quite closely, however (once 
you have found a formula for success, why change alt8 )) 6 


Much of the added volume comes from additional 
examples and case studies, as well as additional exercises. 
The text accompanying the presentation of examples is so 
clear and unambiguous that in many cases I was able to 
anticipate the coded FORTRAN statements prior to looking 
at them. A great deal of the expanded volume comes from 
new material dealing with the interactions of FORTRAN in 
the time sharing environment, although I was disappointed 
not to see an example of a FORTRAN application designed 
specifically for this medium, 


McCracken and his associates are to be congratulated 
for the job they did in assembling and editing the manuscript 
for the second edition. 


The new book is printed on high quality, easy-on- 
the-eyes, non-glare paper and utilizes easily readable 
type. Throughout the book I was able to find only two 
small, insignificant typographical errors and these 
were in areas where you wouldn't ordinarily be looking 
for such errors. The flowcharts used within the book, 
for the most part, follow the ANSI standard. The only 
major criticism I can make of the book is that many of 
the figures and tables, sample programs, etc., appear 
on pages other than those containing the text which 
refers to them. ' 


McCracken's new book has much to offer both the 
neophyte just entering the computing industry and the 
experienced professional who has been around for some 
time. Of greatest importance to beginners is the fact 
that the book will impress him with the fact that FORTRAN 
is a flexible, general purpose problem solving tool and 
is not strictly the toy of the scientist or engineer. 


Book Review 


COMPUTING WITH MINI COMPUTERS 

by Fred Gruenberger and David Babcock 

Melville Publishing Company, Los Angeles, xv + 288 pps. 
$11.95 hardcover 


Reviewed by R. W. Hamming 


The mini computers have arrived. Most of the 
officialdom of computing has long been wrapped up in the 
giant machines, but a few alert people have noticed the 
minis. Fortunately Gruenberger and Babcock have done 
something about the lack of textbooks in computing that 
recognize the mini as a worth while machine on which to 
learn about computing. Their text does not distract 
the beginner from the essence of computing by getting 
a huge mass of hardware plus a vast, ill-documented mess: 
of software systems between the user and computing. 


There are more opinions on how computers should be 
presented than there are people in the field (some of us 
have at least two different opinions). This book adopts 
the down to earth, no "gee-whizz" approach. 


Often the most accurate presentation of a book is 
its table of contents along with the page numbers so that 
the reader can see both the topics covered and the relative 
importance the authors attach to each topic. 


Prologue ab Floating-point arithmetic 113 
Approaching the computer Interpreters 

The sizes of computers Fortran 

Flowcharting More on Fortran 

Binary arithmetic BASIC 

Getting started Some larger problems 

Barriers Program testing 

Assemblers The computing art 

Subroutine packages plus five appendices 


The book is a careful, solid introduction to 
computers and their use in scientific computing (though 
one may wonder how the reader gets to Program Testing, 

page 190, without having to learn it for himself along. 
- the way). The book avoids both business languages and 
business applications, as well as the glamorous area of 
artificial intelligence. 


While the presentation of the material.is based on 
the Varian 620/L, it could easily be adapted to most mini 
machines by supplying the manufacturer's reference 
manuals (which are seldom fit for direct use in teaching 
but are necessary as a reference source for details). 
Thus, because minis are becoming so widespread, this 
book fills a real need. 
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Track 1 (inner) 


Track 2 (middle) 
Track 3 (outer) 


Two markers, A and B, move in the circular pattern 
shown in Figure H. Marker A is about to move three 
sectors forward, to sector 4. Marker B is about to move 
one sector forward, to sector 16. These two moves 
initiate the procedure. 


Subsequent moves are made as follows. Marker A 
is at sector 4, track 1; marker B is at sector 16, 
track 3. A will move on the track specified by B 
(in this case, track 2), which takes A to sector 6, 
track 2. 3B then moves on track 3 (as specified by A) 
from sector 16 to sector 18. For each move, B's 
position dictates the track to be used by A, and A's 
position then dictates the track to be used by B. 


What, then, is the logic of a subroutine which, 
when called, outputs a full move; that is, the next 
position for both A and B? Given that logic, and a 
debugged and tested program for it, it should be possible 
to predict the positions of the markers at any subsequent 
time. To provide a specific goal (and to be able to 
compare independent results), we will seek the positions 
of the two markers (sector and track) at the end of 
1000 full moves. 


The use of the subroutine will make it possible 


to ascertain whether or not A ever overtakes or passes 
B, or vice versa. 


Hewlett-Packard HP-80 


The HP-80 is the "commercial" version of the HP-35. 
In place of the keys that are dear to scientists and 


engineers (logarithms and trigonometric functions), the 


significant functions on the 80 are for compound interest, 
statistical calculations, and a calendar. ‘The machine 
uses floating decimal arithmetic and will expand into 
scientific notation, although there is no mechanism for 
entering numbers in scientific notation directly. . The 
number of decimal places desired in the display (up to 8) 
is programmable. The number of functions available on 
the machine is increased over the HP-35 by use of a shift 
key; that is, several of the keys have two purposes. 


resi 
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The chief utility of the machine lies in its logic 
for compound interest. Of the four factors entering such 
a calculation (the number of time periods; the interest 
rate per time period; the payment per period; and the 
future value), any three may be entered and the machine 
will furnish the fourth. Thus: 


17 at 6%/year for 23 years yields $64.94. 

17 per year at 6%/year for 23 years yields $798.93. 

17 at 6%/year yielding $100, the interest rate is 8.008% 
17 per year for 100 years at 100% yields $2 .1155E31 


Some compound interest calculations can cause the machine 

to cycle for long periods of time. For example, $17 per 

year for 100 years to yield $10000: the interest rate is 
2.95% and the calculation takes 15 seconds. If the future 
value in that calculation is changed to $1,000,000, the 
interest rate is 1431.28% and the calculation takes 30 

seconds. By stepping up the future value and the number 

of time periods, the machine can be made to cycle indefinitely, 
and can only be stopped by shutting off its power. 


For those calculations of compound interest that were 
within range, results were checked against Financial Compound 
Interest and Annuity Tables, Financial Publishing Company, 
1960, 884 pages. 


The calendar buttons will calculate elapsed days 
between two dates within a 200 year range. Thus, entering 
January 1, 1901 and December 31, 1970, the result is 25566 
days. Or, entering January 1, 1901 and 25566, the result 
is December 31, 1970. All the test cases listed on pages 
29 and 33 of Computing: A Second Course (Gruenberger, 1969) 
were in agreement with the results on the HP-80. 


Means and standard deviations are readily calculated; 
the standard deviation formula uses N(N-1) in its denominator 
(4.e., sample deviation). The Trend Line key gives a 
linear regression for data at equally spaced intervals. 


The HP-80 sells for $395. For anyone who deals 
in money calculations, the machine is a marvel. 


LET US ALWAYS STRIVE FOR accurcay @ 
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Canon F-10 


The Canon F-10 ($545) is a 12-digit desk machine with 
a large clear display, weighing 6 pounds and measuring 
10 i/h x 11 1/2 x 3 1/2, AC only. It has floating arith- 
metic (but not scientific notation) plus 0, 1, 2, 3, 4, 5, 
and 6 positions of fixed point. Function keys include 
square and cube root, squaring, common and natural logarithms, 
ex and aX, sine, cosine, tangent, inverse trigonometric 
functions, degrees to radians and radians to degrees, 
rectangular to polar coordinates and the inverse, and 
degrees, minutes, and seconds to decimal degrees. 


All the functions are calculated by iteration. At 
least 10-digit accuracy is claimed for all functions; this 
seems to be correct. The normal tests: 


in x followed by Che (ote produce x 
sin x followed by arcsin x to produce x 
cube root x followed by cube x to produce x 


sin@x + cos’x to produce l. 


all checked to at least 10 significant digits for a wide 
range of x values. Other interesting results were: 


In 2.71828182845 
log 10 
log 100 
log 200000000 
tan71 1 
sin 60 
97, sqrt, sqrt, sqrt, sqrt, 


i. 
dhe 


af 
8.301029996 


5. 
-8660254035 
qrt, Sq, Sq, 8q, Sq, Sq = 
gh 96 99999681 
rectangular V3, 1 to polar = P= 2., 8 = 30.00000001 


TA VE NG 3b bE ae 88 


The trigonometric function keys operate only in degrees; 
sin 36030 = .5, and so on. 


The functionVA was explored at length; see results 
elsewhere in this issue. 


Of the desk-top machines so far reviewed (the Victor 
18-1721 in PC-4 and the Dietzgen ESR-1 in this issue), the 
Canon F-10 is outstanding. It is well designed, reasonably 
priced, and performs exactly as advertised. 
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Desk Colaletor Reyiew 


Dietzgen ESR-1 


The advertising brochure for the Dietzgen ESR-1 
($695) states "every operation, from simple multiplication 
to the most complicated chain calculation, is completed 
in seconds, to 12 digit accuracy." The speed claim is 
excessive; for example, the machine has a cube root key; 
the operation is performed by some interval-halving scheme 
and takes 15 seconds for one cube root. 


Natural logarithm and e* are available; for the 
latter function, the brochure states "Exponent--Automatic- 
ally computes exponent of number displayed." On entering 
the number 3, and then taking the natural logarithm 
followed by the exponential, the result is 2.99999998898, 
which hardly conforms to the claim of 12 digit "accuracy." 
Common logarithms are available, but not the inverse 
function of 10%, The natural logarithm of 2.71828182845 
comes out .996314418; the common logarithm of 10 as 
-999999998. 


The brochure states "Arc Functions--Automatically 
computes arc functions of any number." The only are 
function on the machine is arctangent. 


Again, the key a is described as "Nth power--Raise 
to any power, automatically." The action of the key is 
severely limited; by “any power" is meant integral values 
of n, from 2 to 9. 


Keys are available for hyperbolic sine and cosine; 
for degrees to radians; and for radians to degrees. 
Conversion of degrees, minutes, and seconds to decimal 
degrees is a built-in function. The machine also has 
two storage registers. 


The Dietzgen machine can be switch-programmed for 
floating decimal operation (but not scientific notation) 
or 2, 3, 4, 6, or 8 positions of fixed decimal. The 


- claim is made for "Automatic Round-off--rounds off the 


it 


last digit displayed if five or greater, whatever that 


means. 


The machine is badly designed, inaccurate in its 
calculations, priced outrageously, and touted in fatuous 
terms. A Mark 2 machine is said to be pending, for 
which the Dietzgen company should do its homework. 
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Bees oy the OP Protlen 


In PC3-14 there was given a table of certain 
values for the expression a8, where the values of A 
and B were the square roots of small primes. Some of 
the entries given in that table are reproduced in the 
accompanying chart, recalculated several ways. The 
low order 6 significant digits of the table from PC3 
are compared as calculated on four desk calculators 
and in double precision Fortran. 


2 3 5 a 11 13 
526919 634654 509877 642535 470775 908194 
526918 634654 509877 642535 470774 908194 
526918 634654 509877 642535 470778 908197 
526919 634654 509877 642535 470775 908195 
526919 634655 509879 642537 (470778 908200 


nN 


055549 412396 650954 832911 031865 036964 
055549 412395 650954 832911 031864 036964 
055549 412395 650954 832911 031878 036965 
055545 412377 650953 832909 031855 036962 
055549 412386 650954 832910 031861 036963 


xl 
) 


612109 029437 296751 431002 415774 879364 
612109 029437 296750 431002 415774 879363 
612109 029437 296750 431002 415782 879368 
612110 029434 296753 431008 415774 879363 
612110 029435 296755 431006 415774 879365 


840502 848093 358797 768390 648657 918535 
840502 848093 358797 768389 648656 918535 
840502 848093 358797 768389 648665 918544 
840503 848090 358800 768399 648658 918538 
840503 848092 358804 768401 648660 91854e2 


oe) 
=i 


014455 916446 O62249 357105 474137 968683 
O14454 916446 062248 357105 474137 968682 
O14454 916446 062248 357105 474157 968700 
O1l4454 916441 062248 357104 474137 968678 


bh ~) 
W 2°) 
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Results on an early HP-35 not corrected for the 
defect mentioned on page 6 of PC-1 


Results on an HP-80 
Results on a later model HP-35 
Results on a Canon F-10, reviewed in this issue 


Results in double precision Fortran (to at least 13S) 
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On page 13 of PC-3 there appeared a discussion of 


irrational numbers that lie close to integral values. 
Mr. Herman P. Robinson, of the Lawrence Berkeley Laboratory, 


has pointed out several errors on that page. 


The value of AB, where A and B are irrationals, is 
not necessarily irrational. The simplest counter example 
is A=e, B= ine. 


Further, Mr. Robinson notes that values of irrationals 


can be made arbitrarily close to integral values with ease. 
For example, cos(107%) differs from unity by jo-lo, 


More interesting is the expression (1 +V2_)N, which 
gets closer to an integral value as N increases. The 
table below summarizes some caiculated results. 


—ViI 
A further erratum: the value or igy4 was given 
on page 14 of PC-3 as 131.9997009, but was copied 
incorrectly on page 13. 
The problem of finding values of A and B such that 
‘A lies close to an integer is still open. 


Number of digits Number of 9's 
N to the left of the after the 
decimal point in decimal point. 
(1 A eae 
at 1 .@) 
2 3h ) 
4 2 1 
8 4 3 
16 7 6 
32 13 12 
64 25 24 
128 49 48 
256 98 SNe 
512 196 195 
1024 392 391 


The last value is reproduced on the next page, as 
calculated with 1000-digit precision. é 


The 100 ~ square Tip Prolbler 


PROBLEM ‘] 


In the course of a larger problem, a subroutine is 
needed to traverse the pattern shown. The input to the 
subroutine is a triplet (R,C,N) indicating the Row and 
Column of a cell, and the Number of cells to advance; the 
output is a pair (R,C) indicating the celi at which the 
trip terminates. For any trip that ends outside the 
pattern, the subroutine should furnish the pair (11,10). 
some test cases are: 


10, 
10, 


GR oa aa 


916857 381894481845645133291190504 33384 346182 3250529940377136 
34.4119075342541907081048361452797 15401777 38226872327 89944875 
‘14305635 34335657529154975938905684991897 40400767 8416106442939 
99555847 5666 390154545677 008456626198312155556295338963185799 
267110365897 313197677 9342427 30037 120092270362251009811020913 
61483691467 167530411497 186601524.45055541248594 1439642829 3601 
204069 34959087847 500441905528833.999999999999999999999999999 
999999999999999999999999999999999999999999999999999999999999 
999999999999999999999999999999999999999999999999999999999999 
999999999999999999999999999999999999999999999999999999999999 
999999999999999999999999999999999999999999999999999999999999 
999999999999999999999999999999999999999999999999999999999999 
999999999999999999999999999999999999999999999999999999999999 
999989093178287622853... 
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PC5-14 | 
The logic of the subroutine could be satisfied by 
having a table in storage: 


1 
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and this information could be contained in 101 computer 
words (if the word size is 22 bits or greater), although 
it would probably be more expedient to have it in 303 
words, for ease in manipulating it. 


It would seem that the problem could be attacked 
more logically than by table lookup. A flowchart for 
this subroutine is needed. 


12345678 910 


